Dynamical behavior across the Mott transition of two bands with different bandwidths 
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We investigate tlie role of the bandwidth difference in the Mott metal-insulator transition of a two- 
band Hubbard model in the limit of infinite dimensions, by means of a Gutzwiller variational wave 
function as well as by dynamical mean-field theory. The variational calculation predicts a two-stage 
quenching of the charge degrees of freedom, in which the narrower band undergoes a Mott transition 
before the wider one, both in the presence and in the absence of a Hund's exchange coupling. 
However, this scenario is not fully confirmed by the dynamical mean-field theory calculation, which 
shows that, although the quasiparticle residue of the narrower band is zero within our numerical 
accuracy, low-energy spectral weight still exists inside the Mott-Hubbard gap, concentrated into 
two peaks symmetric around the chemical potential. This spectral weight vanishes only when the 
wider band ceases to conduct too. Although our results are compatible with several scenarios, e.g., 
a narrow gap semiconductor or a semimetal, we argue that the most plausible one is that the two 
peaks coexist with a narrow resonance tied at the chemical potential, with a spectral weight below 
our numerical accuracy. This quasiparticle resonance is expected to vanish when the wider band 
undergoes the Mott transition. 

PACS numbers: 71.30.-|-h, 71.10.Fd, 71.27.+a 



I. INTRODUCTION 



Unlike in single-band models, the Mott metal-insulator 
transition (MIT) in multi-orbital strongly correlated sys- 
tems generically involves other energy scales besides the 
short range Coulomb repulsion U and the bare electron 
bandwidth. They include, for instance, the Coulomb ex- 
change J which produces the Hund's rules, any crystal 
field or Jahn- Teller effect splitting the orbital degener- 
acy, and possibly bandwidth differences between the or- 
bitals. There are many theoretical works making use of 
the so-called Dynamical Mean-Field Theory (DMFT>i 
which analyze the role of the exchange the crystal 
field splitting)^ and the Jahn- Teller effecti^ All these 
analysis suggest that these perturbations, which have the 
common feature of splitting multiplets at fixed charge, 
are amplified near the MIT, leading for instance to an 
appreciable shift of the transition towards lower U^s^ or 
to the appearance of anomalous phases just before the 
MIT.^ This behavior is not surprising, since the more the 
electronic motion is slowed down, i.e., the longer is the 
time electrons stay localized around a site, the larger is 
the chance to get advantage of multiplet-splitting mech- 
anisms. 

On the contrary, the role of different bandwidths for 
nearly degenerate orbitals is less predictable, since the 
Coulomb charge repulsion only cares about the total 
number of electrons at a given site, while it is not con- 
cerned with the orbital they sit in. Recently, this is- 
sue has been addressed in a two-band Hubbard model 
by DMFT, yet leading to controversial results. Liebsch 
has suggested, on the basis of a DMFT calculation us- 



ing quantum Monte Carlo as impurity solver, that both 
with J = and J 7^ the two orbitals undergo a com- 
mon MIT, whatever the difference between their band- 
widths&i This conclusion has been questioned by Koga 
and coworkers)^ again by DMFT, using, however, exact 
diagonalization instead of quantum Monte Carlo. Their 
conclusion is that the MIT is unique only if J = 0. On the 
contrary, when J 7^ 0, there is a first transition at which 
the orbital with smaller bandwidth becomes insulating, 
followed at larger values of the interaction by a second 
transition at which the other orbital ceases to conduct 
as well. This two-stage quenching of the charge degrees 
of freedom has been named Orbital-Selective Mott Tran- 
sition (OSMT) by those authors. Although the coexis- 
tence of localized /-electrons and itinerant d-electrons is 
not unusual in rare-earth compounds, the conclusions of 
Ref.i are a bit surprising in the case of degenerate or- 
bitals. Indeed, the Coulomb exchange-splitting, rather 
than favoring an OSMT, should naively oppose to it, 
since J competes against the angular momentum quench- 
ing due to the different bandwidths. 

In this work, we attempt to clarify this issue by means 
of a variational analysis based on Gutzwiller wave func- 
tions, by standard DMFT calculations as well as by 
an approximate DMFT projective technique. The pa- 
per is organized as follows. In Section ^ we intro- 
duce the two-band model and discuss general properties. 
In Section IIIII we apply a variational technique based 
on Gutzwiller-type trial wave functions to analyze the 
ground state of the Hamiltonian. A full DMFT analysis 
is presented in Section IIVI As a guide to interpret the 
DMFT spectral functions, in Sectional we show the den- 
sity of state obtained by Wilson Numerical Renormaliza- 



2 



tion Group2i of the Anderson impurity model onto which 
the lattice model maps within DMFT. In Section IVTI we 
present an approximate DMFT solution obtained by pro- 
jecting out self-consistently high-energy degrees of free- 
dom, which allow a better low-energy description. Con- 
clusions are drawn in Section IVIII 

II. THE MODEL 

We consider a two-band Hubbard model at half-filling 
described by the Hamiltonian 

2 

w = - E T.*-fLha. + H.c. 

+ ("»"2)'+Wcxch, (1) 

i 

where fj^^^ creates an electron at site i in orbital a = 

1,2 with spin ct, Uia = J2a flaafiMcr t^e occupation 
number at site i in orbital a, and rii = nn + rii i is the 
total occupation number. The explicit expression of the 
Coulomb exchange Ticxch is 

+f E E /m./m.'A2.'A2.+^^-c. (2) 

+ ^ 5^ /'tU(T/?,2f7/i',2cr'/z,lcr' + ^ 

^ 2JJ2 {Tl + Tl), (3) 

i 

where 

Tia — — fi,aa '''ab fi,bu i^) 

a,b a 

are pseudo-spin-1/2 operators, with r" the Pauli ma- 
trices, a = X, J/, z. In what follows, we always take 
< t2 < ^1- Let us start by discussing some general 
properties of this Hamiltonian. 

If C/ 3> ti, the model describes a Mott insulator in 
which two electrons localize on each site. For J > 0, 
the atomic two-electron ground state is the spin triplet, 
followed at energy 2 J by the two degenerate singlets (we 
drop the site index): 

Vl (/i\/2i + /It/li) |0), ^l\ {flA - fhfli) |0): 
and finally at energy 4 J by the singlet: 

\ll{fhfli + fhfn) |0)- 



Hence, the Mott insulator for very large U, specifically 
t1/U ^ J, is effectively a spin-1 Heisenberg model where, 
at any site, each orbital is occupied by one electron, the 
two electrons being bound into a spin-triplet configura- 
tion. Within the OSMT scenario, below some critical 
repulsion, defined in the following as Ui , electrons in or- 
bital 1 start moving, while one electron per site remains 
localized in orbital 2. Only below a lower U2 < Ui, elec- 
trons in orbital 2 dclocalize too. In this particular exam- 
ple with a half-filled shell, the Coulomb exchange does 
not conflict with the OSMT, since Ticxch favors single 
occupancy of each orbital. Yet, one may wonder about 
the role of the pair-hopping term (0) which can trans- 
fer electrons from the delocalized orbital to the localized 
one. 



III. GUTZWILLER VARIATIONAL 
TECHNIQUE 

Let us start by a variational analysis of the ground 
state of the Hamiltonian In particular, we are go- 
ing to use the Gutzwiller variational approach, which is 
one of the simplest way to include electronic correlations 
into a many-body wave function. In what follows, we 
briefly introduce this variational technique for a generic 
multi-orbital model, and later we apply it to our specific 
example. 

Let us consider in general a /c-orbital Hamiltonian that 
contains, besides the hopping term 

k 

^o = - E Y.*-fUha. + H.c., (5) 

{lJ>,<T Q=l 

an on-site interaction of the general form 

^-t = EE'^("'r)^K«.r), (6) 

where Vi{n, F) = |z; n, F) (i; 71, F| is the projector operator 
onto the site-i configuration F with n electrons. If |^'o) is 
the Fermi-sea Slater determinant of the non-interacting 
Hamiltonian, the Gutzwiller wave function j'l') is defined 
through 

l^-) =^g|*o> =n^'G|«'o), (7) 

i 

where the operator Vi g ^icts on site i and is given by 

T'.G = E -^"rPzKF). (8) 

n,r 

The A„r's in (jH)) are variational parameters which mod- 
ify the weights of the on-site configurations |«;n, F) in 
accordance with the interaction © . Following Ref. 0, 
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we assume, without loss of generalityjii that 

(^ol^fcl^o) = (*ol*o> = 1, (9) 
(*0ln'G/L.46.'l*0) = (*oi/Uwi*o> 



N 
2k 



^ab^crcr', 



(10) 



where N is the average occupation per site within the 
Fermi sea, in which all orbitals are equally occupied by 
our choice of the hopping term jSJ. These two conditions 
allow an analytic evaluation of any average value over 
the correlated wave function in infinite dimensions ii£ In 
particular, one can show that Q is satisfied by 



,2 _ Pin.r) 



p(o)(n,r)^ 



where 



P{n,r) = {^\V,{n,T)\^), 



(11) 



(12) 



is the correlated probability distribution of the on-site 
configuration F with n electrons and 



p(«)(n,F) = (vI'o|7',(n,F)|vI/„) 



(13) 



is the uncorrelated one. Therefore, rather than using the 
variational parameters A„r's, one can directly work with 
the correlated probability distribution P(n, F) which sat- 
isfies 

^P(n,F) = l, ^nP(n,F)=iV. 

n,r n,T 

By the above definitions, one can demonstrateiS that the 
average value of the hopping in infinite dimensions re- 
duces to 

= Za{^0\flaJ,,aM, (14) 

where the reduction factor Za can be evaluated through 

^/Z;='^{^0\V^GfU'P^GUaa\'^0)■ (15) 

In addition, the average value of the on-site interaction 
is found to be 

(*|Wint|4')=^^C/(n,F)P(n,F). (16) 

i n,r 



A. Results for J — 

Let us apply this technique to our model |[[J, start- 
ing with the simpler case where J = 0. We define as 
P(0) and P(4) the correlated probabilities of a site occu- 
pied by zero and four electrons, respectively. P(l, -I-) and 
P(l, — ) are instead the probabilities of a singly-occupied 



site with one electron in orbital 2 and 1, respectively. 
Analogously, P(3, -|-) and P(3, — ) are the probabilities of 
a site occupied by three electrons, with one hole in orbital 
1 or 2, respectively. Finally, concerning doubly occupied 
sites, P(2, 1) is the probability of the spin-triplet config- 
uration, P(2, 0) of the inter-orbital spin-singlet, P(2, +) 
and P(2, — ) of the spin singlets with two electrons in or- 
bital 2 and 1, respectively. By particle- hole symmetry, 
we have that 

P(0) = P(4), 
P(l,+) = P(3,-), 
P(l,-) = P(3,+), 
P(2,+) = P(2,-). 

After some simple, but lengthy algebra, we find that the 
hopping energy reduction factors are 



+^/2 v/P(l, -)P(2, +) + V2 ^P{1, +)P(2, 0) 
+V6v/P(l,+)P(2,1), 



% = 2 V2 VP{0)P{1, +) + V2 v/P(l, +)P(2, - 



+V2 v/P(l, +)P(2, -) + V2 v/P(l, -)P(2, 0) 
+\/6V^(l,-)P(2,l). 

Since J = 0, the two-electron configurations with one 
electron in each orbital are equally probable, namely 
P(2, 1) = 3P(2, 0). Therefore, one can use the follow- 
ing parametrization 



^(1,+) 
^(1,-) 
P(2,l) 

P(2,0) 

P(2,+) 



P(l) cos^^ 
P(l) sin^ (f); 

P(2) - cos2 
P(2) - cos^ 



P(2, 



P(2) 



1 



sin^ e. 



The normalization condition now reads 

2P(0) + 2P(1) + P(2)==1, 
and the Z-reduction factors can be written as 



(17) 



2 \/2 sin ^/P{0)P{l) 
+ ^/P{l)P{2) 2 sin sin O + Vs coscj) cos 
2 \/2 cos ^P{Q)P{1) 



2 cos 6 sin ( 



+ /P(T)P(2) 
If we define 



sm cos ( 



i?.c.|*o), 
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the average value per site of the hopping operator in the 
Fermi sea, then the variational energy of the Gutzwiller 
wave function in infinite dimensions is 



E 



~hTZi-t2TZ2 + UPil)+4UP{0). (18) 



Here, the Z^'s are functionals of the probabiUty distri- 
bution P{n), 71 = 0, 1, 2, with the normalization H17() . as 
well as of the two angles 4>,0 € [0, 7r/2]. 

One can proceed analytically a bit further. It is known 
that near a Mott transition and within the Gutzwiller 
wave function approach, one can safely neglect P(0) = 
P(4)ii Within this approximation, P{2) = 1 - 2P(1) 
and 



Zi = ^P{l)P{2) 2 sin(/) sinfl + V8 cos(/) cos6' 



Z2 = VP(1)F(2) 2 cos(/) sin6' + VS sin(^ cos6' 



We denote P(l) = d/2, hence P{2) = \ - d, so that the 
variational energy (|18|) becomes 



where 

am 



2 sin (j) sin + Vs cos (p cos ( 
t2 2 cos cj) sin 6 + VS sin ch cos ( 



The optimal value for d is 

_ Tf{4>,0)-U 

and the variational energy becomes: 



[Tf{^,e)^u] 
8r/(0,0) 



At given cj) and 9, the Mott transition at which both 
orbitals localize occurs when = 0, namely when 

c/,(0,0) = r/(0,0). 

The most stable solution is the one which maximizes /. 
An OSMT corresponds to a situation in which the Mott 
transition occurs with orbital 2 being already strictly 
singly-occupied, namely with (f) ~ 6 ~ 0. This solution 
is an extremum of /. Yet, one has to check whether it 
is also a maximum. We find that this is indeed the case 
whenever 



1 

5 



(19) 



Therefore, within the Gutzwiller variational technique, 
an OSMT can occur even in the absence of Coulomb 
exchange, provided (|19|l is satisfied. 




FIG. 1: Reduction factors Za for the wide band (dashed line) 
and for the narrow band (continuous line) obtained with the 
Gutzwiller wave function for different ratios of the bandwidth 
D2/D1 and J = 0. 



We optimized (|18|l numerically, considering an infinite 
coordination Bcthe lattice at half-filling. The free density 
of state is given by 



Pa(e) 



2nt'i 



(20) 



and, in this case, T = 8/(37r). The half-bandwidth of 
each band is Da = 2ta, and our unit of energy is Di. 
The reduction factors Zi and Z2 arc shown in Fig. ^ 
for different ratios D2/D1. The results indeed confirm 
the analytical calculation H19|) . displaying two distinct 
transitions when D2/D1 < 0.20. 



B. Results for J / 

Let us now move to the more complicated case of J 7^ 
0. The Hund's rule coupling acts only within the two- 
electron configurations and can be written as 

Hj - ^ -|j|i;2,l)(z;2,l| + |j|i;2,0)(i;2,0| 



+ -J(|i;2,-H)(i;2,+| + |i;2,-)(z;2,-|) 

+J (|z;2,+)(*;2,-| + |i;2,-)(i;2,+|), (21) 

using the previous notation for the two-electron config- 
urations. Since it is not diagonal in the representation 
which we have used so far, we are forced to generalize the 
Gutzwiller correlator (jS)) intoii 

r,G ^ V^G + A2± (|2,+)(2,-| + |2,-)(2,+|) . (22) 
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Consequently, 

P(2,+) = A2^pW(2,+) + A^^p(°)(2,-), 
P(2,-) = Ai_p(o)(2,-) + Ai^pW(2,+). 

Since by particle-hole symmetry P(2,+) ~ P(2,— ), as 
well as PW(2,+) = p(")(2,-), then A2+ = X2- ^ X2± 
and hence 



P(2,+) 



2,± p(o)(2,+)- 
In addition, the following matrix element is non zero 
Ai2,± 



so that 

A2,± 



(vI/|2,+)(2,-|vI/) = 2A2,iA2,±P^")(2,+), 

(23) 



2v/P(")(2,+) 



A(2,±) 



+ V^(2,+)- A(2,± 
1 



2/pW(2~^ 



v/P(2,+) + A(2,±) 



^(2,±) 



We notice that |^(2,±)| < P(2,+). The hopping reduc- 
tion factors are consequently modified into 



+V2/P(T~) [/P(2; 



A{2,±) 



V^(2,+)-A(2,±) 



+V2 v/P(f , +)P(2, 0) + V6 v/P(l, +)P(2, 1) 
2V2 VP(0)P(l,-f) 



+V2v/P(f,+) v/^(2,+) + A(2,±) 



+ V^'(2,+)-A(2,±) 



+\/2 v/P(f , -)P(2,0) + \/6 v/P(l,-)P(2, 1). 
The average value of the Hund's coupling is 

Ej = -^JP(2,1) + ^ J [P(2,+) + A(2,±)] 

+ [P(2,0) + P(2,+)-A(2,±)], (24) 

and that of the Hubbard repulsion 

Eu = U [P{l,+) + P{l,-)]+4UP{0), (25) 

thus, the variational energy to be minimized is 

E^-tiTZi-t2TZ2+Ej + Eu. (26) 

By the numerical minimization of the variational energy, 
we find that the critical ratio D2/D1 for an OSMT in- 
creases when J 7^ from the value 0.2 found for J — 0. 



^0^/0^ = 0.40; 




FIG. 2: Reduction factors Za for the wide band (dashed 
line) and the narrow band (continuous line) obtained with the 
Gutzwiller wave function for different ratios of the bandwidth 
D2/D1 and J/U = 0.10. 



u/o 




0.4 0.6 
0^/0^ 



FIG. 3: Phase diagram obtained within the Gutzwiller ap- 
proximation in the absence of the Hund's coupling. Inset: the 
critical ratio D2/D1 below which an OSMT is observed as a 
function of J /U . 



In Fig. 121 we show Zi and Z2 as obtained by numeri- 
cal minimization of (|26|) . for two ratios of D2/D1 and 
J/U = 0.10. Additional calculations allowed to draw 
the phase diagram within the Gutzwiller variational ap- 
proach, see Fig. 13 As is apparent from the inset, the 
introduction of the Hund's coupling increases the value 
of the critical ratio D2/ Di. We also notice that the Mott 
transition at which both bands localize becomes first or- 
der in the presence of J, which, however, might be a 
pathology of the Gutzwiller wave function^ 

In conclusion, we find that the Gutzwiller variational 
technique predicts an OSMT both for J = and J 7^ 0, 
provided D1/D2 is smaller than a critical value which 
increases with J. 



IV. DYNAMICAL MEAN-FIELD THEORY 

To have further insights on the quality of the 
Gutzwiller wave function in infinite dimensions, we have 
performed an extensive DMFT calculation for the same 
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Hamiltonian. DMFT is a non-pcrturbative approach 
that neglects spatial correlations but retains fully time 
correlations and becomes exact for infinite-coordination 
lattices. In this limit, the lattice problem ^ can be 
mapped onto an Anderson impurity model supplemented 
by a self-consistency condition. For simplicity, we con- 
sider an infinite coordination Bethe lattice, as in the 
Gutzwiller variational approach, with a bare density of 
states given by (^0)1 . Again, Da ~ 2ta is the half- 
bandwidth of band a. The Anderson impurity model 
onto which the lattice model maps within DMFT is 

k.a.cT k,a,a 

+ |(n^_2)2-f 2j(T2 + r2^, (27) 

where Tq,, a = x,y,z, are the pseudo-spin operators Q 
for the impurity. The self-consistency condition relates 
the impurity Green's function for orbital a, Ga, to the 
parameters e^a a-nd Vka through 

tlGailiOn) = J2 ■ ■ (28) 

We solve the Anderson impurity model using an exact 
diagonalization methodic at zero temperature. The con- 
tinuous conduction-electron bath is modeled by a finite 
number of parameters eka and Vka {k = 1^ ... ,71^ — 1). In 
our calculations, we considered Us ~ 6 and 4 (not shown). 
The self-consistency (|28|l is implemented through a fit- 
ting procedure along the imaginary axis. To this end, 
we discretize the axis into Matsubara frequencies iujn = 
(2n -I- l)7r//3, where /? is a fictitious temperature that we 
have set to (3 = 500. Moreover, the smallest frequency 
ijJmin has been determined by the smallest pole in the 
continued fraction expansion of the Green's function^ 
Note that this procedure sometimes leads to different 
minimum cutoff frequencies for the two orbitals. In the 
following, we will work in units of Di. 

First, we treat the case in which there is no Hund's 
coupling (i.e., J = 0). The metallic or insulating nature 
of each band is characterized by its quasiparticle residue 
Za'^ = 1 — {dImYia{iio) / diu) |c^^o- In Fig. ^ we show 
Za as a function of the Coulomb repulsion U for different 
bandwidth ratios D2/D1. When the ratio of the band- 
width D2/D1 > 0.20, the quasiparticle weights decrease 
as the Coulomb interaction gets bigger. Even though 
there is a stronger initial reduction for the narrow band, 
the weights eventually vanish for the same critical value 
of Uc/Di ~ 3.6. The situation changes when we further 
decrease the bandwidth ratio (i.e., D2/D1 = 0.15). In 
this case, we find that the weights of the bands vanish 
for different values of U, in agreement with the results of 
the Gutzwiller wave function. Moreover, the critical ra- 
tio D2/D1 = 0.20 that we found earlier seems consistent 
with the DMFT calculation. 

Let us turn now to the model in the presence of a fi- 
nite Hund's coupling, J/U — 0.10, and perform the same 




FIG. 4: Quasiparticle residues Za for the wide band (squares) 
and for the narrow band (circles) obtained by the DMFT 
calculation for different ratios of the bandwidth D2/D1 and 
J = 0. 




FIG. 5: Quasiparticle residues Za for the wide band (squares) 
and the narrow band (circles) obtained by the DMFT cal- 
culation for different ratios of the bandwidth D2/D1 and 
J/U = 0.10. 



calculations for the following ratios of the bandwidths: 
D2/D1 = 0.80 and 0.40. The results are shown in Fig.|Sl 
We still find evidences for an OSMT, this time, however, 
below a larger ratio of the bandwidths. Further calcu- 
lations show that the critical ratio of the bandwidth for 
J/U = 0.10 is D2/D1 - 0.60. 

If we were to confine our analysis to the behavior of 
the quasiparticle residues Z^'s, we should conclude that, 
both in the absence and in the presence of a Hund's 
coupling, the OSMT scenario does occur for sufficiently 
small ratio of the bandwidth, in qualitative and also 
quantitative agreement with the variational results of the 
Gutzwiller wave function. The only difference may lie in 
the order of the transition. Even if it is difficult to set- 
tle precisely the order of the transition within DMFT, 
our results seem to point towards a second-order phase 
transition at the MIT Ui with finite J, contrary to the 
first-order transition predicted by the Gutzwiller wave 
function. 

A deeper insight into the above scenario can be gained 
by analyzing the spectral properties of the more corre- 
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FIG. 6: Density of states for the wide band (left panels) 
and the narrow band (right panels) as obtained by DMFT for 
different values of the Coulomb interaction. The ratio of the 
bandwidth is D2/D\_ = 0.15 and J = 0. 
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FIG. 7: Upper panels: low-energy part of the density of 
states of the narrow band as obtained by DMFT for different 
values of J/Di. The ratio of the bandwidths is D2/-D1 = 0.15, 
and the Coulomb repulsion is set to U / D\ =2. Lower panel; 
position of the low-energy peaks as a function of J/Di . 



arc also present at J = 0, and move linearly away from 
the chemical potential when J ^ 0, roughly as 2J, see 
Fig. Their total spectral weight scales approximately 
like the quasiparticle residue of the wider band, Zi, both 
vanishing at the second MIT, Ui. In addition, if J = 0, 
the distance between the peaks also scales like Zi, see 
Fig. El If we, reasonably, assume that these two peaks 
mimic two resonances, one below and the other above the 
chemical potential, it becomes much less obvious what 
might be the actual value of the DOS right at the chem- 
ical potential if we were not constrained to a small num- 
ber of levels. Moreover, even if the DOS were strictly 
zero at the chemical potential, still we should determine 
whether band 2 behaves like a small-gap semiconductor 
or a semimetal for U2 < U < Ui. In other words, the en- 
ergy discretization inherent in the exact diagonalization 
technique might play a more critical role in our case than 
in the simplest single-band Hubbard model. 

Therefore, although the numerical evidences we have 
presented so far point in favor of the existence of an 
OSMT with zero or finite J below a critical bandwidth 
ratio, there are several aspects which still need to be 
clarified. We will consider a deeper investigation of such 
aspects in the following sections. 



V. SINGLE IMPURITY SPECTRAL 
PROPERTIES 

The first issue we want to address concerns the origin of 
the two peaks in the orbital 2 spectral function inside the 
Mott-Hubbard gap. Although the self-consistency condi- 
tion (|28|) of the effective Anderson impurity model H27|) 
plays a very crucial role, for instance it determines a criti- 
cal value of U above which the Kondo effect docs not take 
place anymore, yet useful information can be obtained by 
studying 1)27(1 without imposing (|28|l . which is what we 
are going to do in this section by means of the Wilson 
Numerical Renormalization Group (NRG)iS 

The Anderson impurity model (|27|l is controlled by sev- 
eral energy scales, the Hubbard U, the Himd's coupling 
J and the so-called hybridization widths 



lated band, and not just its quasiparticle residue. In- 
deed, such an inspection leads to a less clear-cut pic- 
ture, revealing features which are not captured by the 
Gutzwillcr wave function. In Fig. wc show the den- 
sity of states (DOS) of both orbitals for various Hubbard 
t/'s and D2/D1 = 0.15. We notice that, although the 
DOS of the narrow band right at the chemical potential 
becomes zero within our numerical accuracy above U2, 
there is still low-energy spectral weight inside the Mott- 
Hubbard gap. Due to our discretization procedure, this 
weight is concentrated in two peaks located symmetri- 
cally with respect to the chemical potential. These peaks 



For simplicity, we will assume that the two conduction 
baths are degenerate with half-bandwidth D, which will 
be our unit of energy. In Fig. |S1 we show the impurity 
spectral function of the orbital 2, ^2(0;), as obtained by 
NRG for J 0, [/ = 2D, Fi ^ D/{2tt), and for several 
values of r2/ri < 1. Since we do not impose any self- 
consistency, the DOS for any r2 7^ shows a Kondo 
resonance at the chemical potential, which narrows as r2 
decreases. In addition, there are two more peaks which 
move slightly away from the chemical potential as r2 is 
reduced. These peaks actually resemble those we find 
in the DMFT calculation. Indeed, they move linearly as 
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FIG. 8: Spectral function A2{uj) of the narrow band for J = 
0, U/D = 2, Ti = D/2TV, and Fi/Ti = 1,1,1,1,0 (from 
top to bottom). Inset: the same spectral function on a wider 
scale. 



1.25 



0.75 




0.25- 



0.05 



(O/D 



FIG. 9: Spectral function A2{iL>) of the narrow band for 
J/D = 0, 0.004, 0.008, 0.012 (from top to bottom), U/D = 2, 
Fi — D/2-R, and F2/ri = 1/8. Inset: the same spectral func- 
tion on a wider scale. 



we switch on J, see Fig. O just like we observe within 
DMFT. 

The origin of these peaks is easy to identify when 
r2 = 0. When the orbital 2 is not hybridized with its 
bath, its occupation number n2 is a conserved quantity. 
The ground state is expected to belong to the subspace 
with 77.2 = 1, because, in this case, (ni) = 1 and the 
Kondo-screening energy gain is maximum. This state is 
twofold degenerate reflecting the free spin- 1/2 of the elec- 
tron localized in orbital 2. The energy gap to the lowest 
energy states for n2 = 0,2, (rii) = 2,0, respectively, is 
therefore of the order of the Kondo temperature of or- 
bital 1, Tki- The impurity orbital-2 DOS is analogous 



to the core-hole spectral function in X-ray absorption, so 
it should start above a finite threshold proportional to 
Tki- In other words, the DOS has a small but finite gap 
of order Tki, similar to what we observe within DMFT. 

However, as soon as r2 is non zero, this gap is filled 
and, in addition, a Kondo-resonance appears. Even if 
we move the peaks away from the chemical potential by 
increasing J, see Fig. the region between them and 
the narrow (practically invisible in the figure) Kondo- 
resonance is still covered by spectral weight. In the light 
of this dynamical behavior, it is not at all obvious what 
the self-consistency requirement (j28|l may lead to when 
band-1 is still conducting. In other words, either a true 
narrow gap, as if r2 = 0, or a pseudo-gap with a power- 
law vanishing DOS, or two-peaks plus the narrow res- 
onance are equally compatible with the self-consistency 
condition. However, the event in which most of the spec- 
tral weight is concentrated in the two symmetric peaks, 
leaving only negligible weight within the narrow reso- 
nance, is extremely hard to identify with a limited num- 
ber of levels. As an attempt to discriminate among the 
aforementioned possible scenarios, in the following sec- 
tion, we implement a projective self-consistency tech- 
nique which allows a more detailed low-energy descrip- 
tion within DMFT. 



VI. PROJECTIVE SELF-CONSISTENT 
TECHNIQUE 



A remarkable feature uncovered by DMFT nearby a 
MIT is the clear separation of energy scales between well 
preformed high-energy Hubbard bands and lingering low- 
energy itinerant quasiparticles. It has been shown^ that 
this partition of energy scales allows to reformulate the 
problem into a new one, in which the high-energy part is 
projected out. Essentially, the original Anderson impu- 
rity model, which involves both high-energy side-bands 
and low-energy quasiparticles, is reduced to a Kondo-like 
model which can be attacked more easily by a numerical 
procedure. In this section, we apply a projective tech- 
nique to our model, which is similar to Ref. Il3l with the 
only difference that the resulting effective problem is still 
an Anderson impurity model with rescalcd parameters. 

As we showed, the occurrence of an OSMT does not 
seem to require a finite exchange but rather a sufhciently 
small bandwidth ratio. Therefore, we prefer to present 
the projective technique in the simpler case where J = 0. 
Following Ref. 0, we start by rewriting the Anderson 
impurity model 1)27(1 explicitly separating low- {L) and 
high- {H) energy scales 



fiAM = ri + At + /T 



(29) 
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where 

+ ^"Aa^Cka. + yiiifl^ka. + ^.c), (30) 

k.a.cr k.a.cr 

describes the impurity coupled to the high-energy levels, 

-H"- = Y."^ 4acL.C,,,. (31) 

k,a,cr 

is the low-energy bath Hamiltonian, and finally 

= ^fca iflc,,. + h.c.), (32) 

k,a,a 

mixes low- and high-energy sectors. The impurity 
Green's function, also written as sum of a low- and a 
high-energy part, Ga{i(jj) = G^{iu}) + G^{iuo), should 
satisfy the self-consistency requirement (|28f) . If we as- 
sume that the low-energy spectral weight is Wa <C 1, 
the self-consistency condition for the integrated low- and 
high-energy spectral functions, p^{() and p^(e), respec- 
tively, implies the following sum-rules 

(^^-)' = ^'W^- (33) 

k 

Yl" i^"af = tl{l-Wa), (34) 

k 

showing that the impurity is strongly hybridized with the 
high-energy levels and very weakly with the low-energy 
ones. Let us for the moment neglect the coupling to the 
latter. The ground state of (|30|l is the adiabatic evolution 
of the states in which all negative-energy bath levels are 
doubly occupied and two electrons sit on the impurity, 
giving rise to a six-fold degenerate ground state. Other 
states with the same number of electrons lie above the 
ground state at least by an energy U . The lowest-energy 
states with one more (less) electron are more degenerate, 
since they emerge adiabatically from the states obtained 
by adding (removing) an electron either in the impurity 
levels or in the positive (negative)-energy baths. This 
large degeneracy is, however, split linearly by V^^, which 
implies the broadening of the Hubbard bands around 
their centers of gravity ±11/2. The main effect of the 
mixing term H32|) is to provide a Kondo exchange cou- 
pling between the six-fold degenerate ground state of 
and the low-energy baths, which can be obtained 
by degenerate second-order perturbation theory in Ti™ 
or, more formally, by a Schrieffer- Wolff canonical trans- 
formationpiSiii Once the effective Kondo model is ob- 
tained, we could for instance follow Ref. lTsL namely solve 
that model and impose the self-consistency condition to 
the impurity Green's function, calculated through the 
Schrieffer- Wolff canonically transformed faa- To be con- 
sistent, one should in principle expand the transformed 



facr up to second order in V^/U and impose the self- 
consistency requirement in the whole energy range, in- 
cluding low and high energies. In practice, even if the 
self-consistency is imposed only to the low-energy spec- 
trum, one still gets a faithful description of the critical 
behavior near the MIT>i^ An equivalent procedure, that 
we have instead decided to follow, is to identify a new 
two-orbital Anderson impurity model, coupled only to 
the low-energy levels, which maps to the same Kondo 
model and next impose the self-consistency only to the 
low-energy part of the impurity Green's function: 



k 




Regarding the high-energy part of the self-consistency, 
since we always model the high-energy levels with just 
four levels at energies = ±C//2, we need to impose an 
additional requirement besides 1)35(1 . which, through H34|l . 
is simply 

Vli = ta y^^, (36) 

where Wa is the low-energy spectral weight obtained self- 
consistently from 1)35(1 . The advantage of the projective 
method is that we can now model the low-energy spec- 
trum with more levels, the cost being the additional self- 
consistency condition 1(36(1 . 

When we apply this projective technique to our two- 
orbital model with J = and ti >t2, we end up with an 
effective Anderson impurity model 

Heff = Y^n,- 1)2 + ^in.2 ir + ^("1 - 1)("2 - 1) 

+ E^ 'Lcl^^kaa + E"^ ytiflc,,, + h.c), (37) 
k,a,a k,a,a 

where Ui , and U2 arc found from the solution of the high- 
energy problem Moreover, U12 = {Ui + U2)/2, 
which assures the six-fold degeneracy of the isolated 
doubly-occupied impurity and we have that Ui < U2 with 
U2 — Ui ^ ti — t2- In other words, the high-energy lev- 
els provide a partial screening of the Hubbard repulsion, 
more efficient within orbital 1, which is more hybridized 
with the bath. Therefore, the difference of bandwidths 
acquires in our projective method a quite transparent 
role: while the bare Coulomb repulsion does not care 
about the orbitals in which electrons sit, this indifference 
is lost once the high-energy screening is taken into ac- 
count. In Fig. IIUI we compare the imaginary part of the 
Green's functions in Matsubara frequencies ll}„ as func- 
tion of LUn as obtained by full DMFT or using the above 
projective self-consistent technique (PSCT), at J = 0. 
Note that within the PSCT, we model the low-energy 
conduction bath through five discrete levels. The agree- 
ment is satisfying, and the additional levels clearly allow 
for a more accurate description of the low-energy Green's 
function. 
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FIG. 10: Imaginary part of the Green's function for D2/ D\ — 
0.18, U/Di = 2.00, and J = as obtained by DMFT and with 
the PSCT. 
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FIG. 12: Low-energy part of the density of states of the 
wide band (left panels) and the narrow band (right panels) 
obtained within the PSCT for different U/Di. The ratio of 
the bandwidth is D2/D1 = 0.18 and J = 0. 




U/Dj 

FIG. 11: Quasiparticle residues Za (triangles) and low-energy 
spectral weights Wa (circles) as obtained by the PSCT with 
J = 0. Dashed lines refer to the wide band, and continuous 
lines to the narrow one. 



In Fig. ^] wc plot the PSCT values of the quasipar- 
ticle residues Zi and Z2 as function of U at J = for 
D2/D1 = 0.18, as well as of the full spectral weights, 
Wi and W2, inside the Mott Hubbard gap. In agree- 
ment with standard DMFT, we find a region where Z2 
is zero within our numerical accuracy, while Zi is still fi- 
nite. Yet, the total spectral weights are both non zero. In 
Fig. El we draw the low-energy DOS for the two bands 
and various C/'s. 



VII. DISCUSSION AND CONCLUSIONS 

In this work, we have studied by several techniques 
the properties of the Mott transition in an infinite- 
dimensional Hubbard model with two bands having the 
same center of gravity but different bandwidths, both in 
the presence and in the absence of the Hund's exchange 
splitting J. We have shown that a variational calculation 
based on a Gutzwiller wave function predicts that the two 
bands may undergo different metal-insulator transitions 
both for J = and J 7^ 0: by increasing U, the narrower 
band ceases to conduct before the wider one. The nec- 
essary condition for this orbital-selective Mott transition 
is that the bandwidth ratio is lower than a critical value 
which increases with J, being 0.2 when J = 0, see Fig.O 
This result would contradict both Liebschji^ who claims 
that there is only a single Mott transition whatever is J 
and the bandwidth ratio, as well as Koga and coworkers;^ 
who instead suggest that the Mott transition is unique if 
J = but splits into two distinct ones as soon as J 7^ 
for any value of the bandwidth ratio. The behavior of the 
quasiparticle residues as obtained by DMFT using exact 
diagonalization confirms, even quantitatively, the varia- 
tional results, showing that the residue of the narrower 
band may vanish before the one of the wider band if the 
bandwidth ratio is sufficiently small, both for J = and 
J 7^ 0. Actually, an OSMT which occurs both in the 
absence and in the presence of an exchange splitting is 
somehow more conceivable, because the role of J is model 
dependent. We notice that, in more general situations 
where the number of orbitals is greater than two and dif- 
ferent from the number of electrons, as for instance in the 
case of t2g orbitals occupied by two or four electrons on 
average, the Coulomb exchange would instead compete 
against the angular momentum quenching which occurs 
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in the OSMT scenario. Therefore, we suspect that the 
role of the Coulomb exchange might actually depend on 
the specific model. 

However, a closer inspection to the low-energy spec- 
tral properties of the narrower band in the region where 
it is apparently insulating while the wider band still con- 
ducts poses doubts to the above simple scenario. The 
reason is that, in spite of a quasiparticlc residue which is 
zero within our numerical accuracy, the narrower band 
has spectral weight inside the Mott-Hubbard gap, which 
scales like the quasiparticlc residue of the wider band. In 
other words, the charge fluctuations which still occur in 
the wider band are transferred into the narrower one, as 
somehow predictable. This low-energy spectral weight is 
concentrated in two peaks symmetrically located around 
the chemical potential. Roughly speaking, the distance of 
each peak from the chemical potential is 2 J plus a quan- 
tity of the order of the quasiparticlc resonance width of 
the wider band. Due to our limited numerical resolution, 
we can not establish rigorously whether these two peaks 
(a) signal a narrow-gap semiconducting behavior, (b) sig- 
nal a semimetallic behavior, with a power-law vanishing 
density of states, (c) or coexist with an extremely nar- 
row resonance at the chemical potential, with a spectral 
weight well below our numerical accuracy, just like the 
single impurity does. Although the elements at our dis- 
posal do not definitely allow to discriminate among these 
three scenarios, yet one can recognize that some of them 
are more plausible than the others. 

The first possibility (a) of a narrow-gap semiconduc- 
tor seems very unlikely. Indeed, in this case, the gap 
between the two low-energy peaks would open large and 
then diminish as the quasiparticlc resonance width of the 
wider band, by further increasing the repulsion U. There- 
fore, the insulating character of the narrower band would 
weaken by increasing U, which seems a bit odd. 

Let us consider instead the scenario (b) of a semimctal. 
If taken literally, it would imply a vanishingly small lo- 
cal magnetic susceptibility, while we actually find a very 



large one, much larger than the local susceptibility of the 
wider band. However, a semimetallic behavior would im- 
ply, in our particle-hole symmetric case, a breakdown of 
Fermi liquid theoryii^ Therefore, a power-law vanishing 
single-particle DOS might not necessarily confiict with 
almost free-spin excitations in a scenario in which Fermi 
liquid theory breaks down and for instance spin-charge 
separation emerges. Although it might represent a quite 
interesting circumstance, yet we could not find any phys- 
ical arguments justifying such a non-Fermi liquid behav- 
ior. Therefore, we are tempted to discard it in favor of 
the more conservative scenario (c) in which the two peaks 
coexist with a narrow resonance which remains tied at 
the chemical potential, its spectral weight being smaller 
than our numerical accuracy. This resonance should dis- 
appear right at the same U where the wider band ceases 
to conduct. 

In any case, the rich structure of the low-energy spec- 
trum, revealed by our study, highlights the subtlety of 
the present problem and could explain the discrepan- 
cies between previous studies. Further investigation is 
needed to confirm our speculation that the narrow reso- 
nance peak displayed by the single impurity survives the 
DMFT self-consistency. 
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